Kinetic Study and Reaction Mechanism of the Gas-Phase Thermolysis Reaction of Methyl Derivatives of 1,2,4,5-Tetroxane

Tetroxane derivatives are interesting drugs for antileishmaniasis and antimalaric treatments. The gas-phase thermal decomposition of 3,6,-dimethyl-1,2,4,5-tetroxane (DMT) and 3,3,6,6,-tetramethyl-1,2,4,5-tetroxane (acetone diperoxide (ACDP)) was studied at 493–543 K by direct gas chromatography by means of a flow reactor. The reaction is produced in the injector chamber at different temperatures. The resulting kinetics Arrhenius equations were calculated for both tetroxanes. Including the parent compound of the series 1,2,4,5-tetroxane (formaldehyde diperoxide (FDP)), the activation energy and frequency factors decrease linearly with the number of methyl groups. The reaction mechanisms of ACDP and 3,6,6-trimethyl-1,2,4,5-tetroxane (TMT) decomposition have been studied by means of the DFT method with the BHANDHLYP functional. Our calculations confirm that the concerted mechanism should be discarded and that only the stepwise mechanism occurs. The critical points of the singlet and triplet state potential energy surfaces (S- and T-PES) of the thermolysis reaction of both compounds have been determined. The calculated activation energies of the different steps vary linearly with the number of methyl groups of the methyl-tetroxanes series. The mechanism for the S-PES leads to a diradical O···O open structure, which leads to a C···O dissociation in the second step and the production of the first acetaldehyde/acetone molecule. This last one yields a second C···O dissociation, producing O2 and another acetone/acetaldehyde molecule. The O2 molecule is in the singlet state. A quasi-parallel mechanism for the T-PES from the open diradical to products is also found. Most of the critical points of both PES are linear with the number of methyl groups. Reaction in the triplet state is much more exothermic than the singlet state mechanism. Transitions from the singlet ground state, S0 and low-lying singlet states S1–3, to the low-lying triplet excited states, T1–4, (chemical excitation) in the family of methyl tetroxanes are also studied at the CASSCF/CASPT2 level. Two possible mechanisms are possible here: (i) from S0 to T3 by strong spin orbit coupling (SOC) and subsequent fast internal conversion to the excited T1 state and (ii) from S0 to S2 from internal conversion and subsequent S2 to T1 by SOC. From these experimental and theoretical results, the additivity effect of the methyl groups in the thermolysis reaction of the methyl tetroxane derivatives is clearly highlighted. This information will have a great impact for controlling these processes in the laboratory and chemical industries.


Introduction
Organic peroxides are a series of compounds with many different applications, due to their therapeutic properties and explosive character.Some of them are used to fight malaria in subtropical regions, and as antibiotic; some others therapeutical properties are the intervention in biological metabolites, cancer illness, cellular aging, and herbicide use [1][2][3][4][5].
Furthermore, commercial vinyl polymers-such as polystyrene, polyethylene, and polyvinyl chloride-are produced industrially through free radical polymerization, where organic peroxides are commonly used as free radical initiators.The thermal decomposition of organic diperoxides through the homolytic cleavage of the peroxide bond generates a biradical intermediate, which would initiate the polymerization of vinyl monomers.This biradical species would be incorporated into the polymer structure, giving polymeric species that contain peroxide groups in their structures, which are then decomposed in the next processes.From the structure of the organic peroxide, a biradical is generated, which influences on the properties of the polymer.Therefore, it is important to clarify the decomposition reaction mechanism of the di-and triperoxide organic compounds, which can be considered multifunctional initiators due to the presence of more than one peroxide group [1].However, the knowledge of the behavior and reaction mechanism of these compounds is limited.Hence, a further understanding of these processes would allow the control of many industrial processes with high impact in economic costs, and environmental effects.
The tetroxane-derivative antimalaric drugs are compounds with large molecular weight, low toxicity and no genetic toxicity.It is well known that the diperoxide compounds are more suitable as drugs than the equivalent ozonides and trioxane compounds [2].Therefore, the diperoxide compounds are considered as new anti-malaric drugs and can be a good choice in front of artemisinin derivative drugs [2][3][4].A deeper knowledge of the reactivity of these compounds close to receptor environment would help in the design of optimal antimalaric drugs.
In previous papers, the experimental kinetic and theoretical mechanism of the thermolysis reaction at gas phase of the parent compound of the series (1,2,4,5-tetroxane (formaldehyde diperoxide, FDP)) was accomplished [7].Additionally, the theoretical mechanism of the thermolysis reaction of 3-methyl-1,2,4,5-tetroxane (MFDP) [8] and 3,6-dimethyl-1,2,4,5-tetroxanes (DMT) [9] was also determined.Previously, a concerted mechanism of the thermolysis reaction was proposed [10], but from quantum mechanics calculations, a stepwise mechanism was also determined.In these previous works [7][8][9], a common stepwise mechanism for the thermolysis reaction was found with lower activation energy than the previously proposed concerted mechanism.In the stepwise mechanism, the first step brings the reactant to a diradical open (o) intermediate in its singlet state (S 0 ), and, afterwards two step-wise mechanisms follows: (i) a decomposition step by step from the o structure to two aldehydes molecules plus an oxygen molecule in its excited singlet state (S 1 ); and (ii), a transition from the S-o diradical structure to an o structure in its triplet state (T 1 ), and a subsequent triplet stepwise mechanism, producing two molecules of aldehydes plus one oxygen molecule in its ground triplet state (T 1 ).This singlet-to-triplet-state non-adiabatic transition is deeply studied in this work.In general, the rate-limiting step was the second step.The thermal decomposition mechanism in dioxetane was also found to proceed in a stepwise manner and with a non-negligible nonadiabatic transition from the singlet to triplet state [11].
The methyl group substitution effect was analyzed previously [9], showing that the activation energies (Ea), reaction energies (Er), and reaction coordinate lengths of the different steps of the stepwise mechanism of the reaction behave linearly with the number of methyl groups.
Because of the peroxide bonds and the cyclic structure, these compounds become highly explosive.In general, they are compounds of an unstable and explosive nature due to the peroxide groups, which are intrinsically unstable and highly reactive [5], so the synthesis and manipulation of these materials require strict care to avoid accidents.The multiple incidents recorded in recent decades involving this type of species have raised the need for a constant and exhaustive physicochemical study.Methyl derivatives are especially dangerous because of their high explosivity, and their synthesis and manipulation should be performed in very small amounts.Additionally, the products formed during the thermolysis reaction of the tetroxane derivatives augment three times the number of moles with respect to the reactants, contributing also to their high-explosive character.The successive methylation of the tetroxanes yields different reaction heats.These features plus the surplus incentive of knowing the kinetic and mechanism of thermolysis reaction of the highest and lowest methyl-substituted members of the series and its comparison with the parent compound increase the interest of this work.The knowledge of these complex non-adiabatic processes will facilitate the disposition of control techniques for further synthesis and manufacturing scale processes.This will have a great economic and environmental impact in the chemical industry.
Herein, the following goals are pursued: (i) to carry out an experimental kinetic characterization of DMT and ACDP of the thermolysis reaction in the gas phase; (ii) to determine the theoretical thermolysis reaction mechanism of trimethyltetroxane (TMT) and ACDP; and (iii) to theoretically revisit tetroxane and its methyl derivatives in order to study the non-adiabatic transition occurring in the diradical open structure (chemiexcitation) by multiconfigurational quantum chemistry.The obtained data will be studied in order to know the methyl substitution effect in the gas phase thermolysis reaction of the methyl-tetroxanes.

Analytical Results
After the introduction of ACDP [10] into the injector chamber [12], only three compounds were detected: solvent (retention time = 1.3 min), acetone (2.84 min), and the remaining non-reacted ACDP (8.7 min).As a consequence of the very short residence time (40 s) of ACDP in the injector chamber, acetone was the only product detected, and no other products coming from possible secondary reactions were found.
The introduction of DMT into the injector chamber showed the detection of three compounds: acetaldehyde (1.4 min), the solvent (3.2 min), and the remaining non-reacted DMT (7.2 min).Acetaldehyde was the only product detected, a consequence of the very short residence time (40 s) of DMT in the injector chamber.

Kinetic Results
According to Arrhenius's equation, the Ea was calculated from the slope of least-square fitting of ln[ln(C 0 /C)] as a function of 1/T and the intersect of this fitting yields the ln(A.t).Considering the retention time, which is constant in all measures, A can be obtained.
The rate constant of reaction, k exp , for DMT and ACDP at different temperatures (   Therefore, the Ea values are 26.5 ± 1.0 and 22.7 ± 1.0 kcal/mol for DMT and ACDP, respectively, and 1.4 × 10 11 and 5.8 × 10 9 s −1 for the frequency factors of DMT and ACDP, respectively.In a previous work, Ea for ACDP turned out to be 39 ± 2.5 kcal/mol [13], coming from a very different technique, where many reactors (Pyrex vessels) were introduced at a constant temperature (at different temperatures), and the reaction was stopped at different times (up to 50 h) without controlling that all species were in the gas phase.Other cyclic peroxides, such as dioxetane (Ea = 22.7 kcal/mol) [11] dimeric diphenyl tetroxane [12] and FDP [7] showed Ea = 21.6 and 29.3 kcal/mol, respectively.From these results, it is possible to observe a substitution effect on the Ea of this series of compounds.The methyl substitution from FDP to dimethyl tetroxane (DMT) produced a decrease of Ea, where the gas phase of DMT is made of a mix of the different position isomers [9].It is possible to infer that the methyl substitution would go on affecting the Ea of these derivatives.Ea values for these three derivatives (FDP, DMT, and ACDP) are fitted as a function of the number of methyl groups (n), resulting in the following linear function: Ea = (29.28± 0.04) − (1.37 ± 0.01)n; R 2 = 0.99 (3) This indicates an average decrease of 1.4 kcal/molCH 3 for Ea.The average value of the frequency factor of FDP was 5.2 × 10 13 s −1 .Our frequency factor shows a much smaller value, having a parallelism with the methyl substitution of the dioxetanes [14,15], which were in the range of 2 × 10 12 -1.3× 10 13 s −1 .Nonetheless, a new linear relationship is also found with the number of methyl group: Therefore, we have to emphasize the additivity of the methyl groups in the thermolysis reaction of the methyl tetroxane derivatives, which could bring further to think that the additivity of the methyl groups in the properties of these compounds would be strongly feasible.The hyperconjugative effect of methyl groups increases the delocalization of the electronic cloud stabilizing the transition states and decreasing the activation energy.
Previous experimental and theoretical works [7][8][9] indicated a stepwise mechanism.In this case, we are going to undertake a computational study in order to confirm and widen the previous proposed mechanisms.

Thermolysis Reaction Mechanism
The possible reaction mechanisms are given in Scheme 1 and analyzed with our calculations, in which we can describe a concerted mechanism (S-I) and a stepwise mechanism

Thermolysis Reaction Mechanism
The possible reaction mechanisms are given in Scheme 1 and analyzed with our calculations, in which we can describe a concerted mechanism (S-I) and a stepwise mechanism (S-II) via a diradical open

Concerted Mechanism
The tri/tetra-methyl-substituted-tetroxane reactants are in a chair conformation (S-c).The H atoms of the methyl groups in geminal position show a symmetric eclipsed conformation two to two with respect to the average plane of the molecule, in which all carbons are included (Figure 1a).One H atom of each methyl group is in this plane, and two H

Concerted Mechanism
The tri/tetra-methyl-substituted-tetroxane reactants are in a chair conformation (Sc).The H atoms of the methyl groups in geminal position show a symmetric eclipsed conformation two to two with respect to the average plane of the molecule, in which all carbons are included (Figure 1a).One H atom of each methyl group is in this plane, and two H atoms out of the plane of two geminal methyl groups are facing each other.The concerted mechanism has been calculated as one single step, where the transition state (TS) is called S-TS X (Table 2, Figure 1a), with a transition vector (TV) of opening and closing the peroxide O-O bonds and opening C•••O bonds to form the O 2 molecule and two molecules of acetaldehyde/acetone via the opening/closing of the C-O bonds.This TS and TV are depicted in Figure 1a.The Ea values are 66.3, 68.8, and 65.9 kcal/mol for TMT-axial (TMTax), TMT-equatorial (TMTeq), and ACDP, respectively, being lower than the less-substituted compounds of the series [7][8][9].The TMTax has lower Ea than TMTeq probably due to the anomeric effect that stabilizes the axial isomer TS.
that the Ea decreases linearly with the increasing number of methyl groups, with a slope of −2.04 kcal/molCH3 (Figure 1b; in this figure, average values of axial and equatorial isomers of TMT are also considered).This slope is larger than that of the experimental Ea as a function of the number of methyl group (Equation ( 3)).This fitting shows that the FDP Ea for a concerted mechanism would be 73.8 kcal/mol, being considerably higher than the experimental value (29.3 kcal/mol) [7].Additionally, the value of Ea = 65.9 kcal/mol for ACDP is the lowest of the series, which is much higher than our experimental value but is still high enough to avoid the concerted mechanism in front of the stepwise mechanism (vide supra).Taking into account all methyl derivatives, a correlation of the average values of Ea as a function of the number of methyl groups has been performed (Figure 1b).We observe that the Ea decreases linearly with the increasing number of methyl groups, with a slope of −2.04 kcal/molCH 3 (Figure 1b; in this figure, average values of axial and equatorial isomers of TMT are also considered).This slope is larger than that of the experimental Ea as a function of the number of methyl group (Equation ( 3)).This fitting shows that the FDP Ea for a concerted mechanism would be 73.8 kcal/mol, being considerably higher than the experimental value (29.3 kcal/mol) [7].Additionally, the value of Ea = 65.9 kcal/mol for ACDP is the lowest of the series, which is much higher than our experimental value but is still high enough to avoid the concerted mechanism in front of the stepwise mechanism (vide supra).

Structure
The reaction product (S-P) (Branch S-I, Scheme 1) shows an O 2 molecule in the singlet state surrounded by hydrogen bonds of two acetaldehyde/acetone molecules depending on if the reactant is TMT or ACDP.Reaction energy (Er) is exothermic, being −31.1 and −28.0 kcal/mol for TMT axial and equatorial isomers, respectively, and −32.8 kcal/mol for ACDP.The series follows a linear equation as a function of the number of methyl groups: Er = −17.1 − 4.0n (R 2 = 0.98) (FDP is not included in this equation).

Stepwise Mechanism in the Singlet State
The stepwise mechanism is depicted in Scheme 1 (Branch S-II).The mechanism has three steps: (i) the production of the open diradical structure (S-o), (ii) the formation of acetone/acetaldehyde (S-b ′ ) and oxide-peroxide intermediates (S-b), and (iii) the generation of the end products (S-p), where the oxygen molecule is in the singlet state.Step (i) is common to all methyl derivatives of tetroxane.In Step (ii), only one possible reaction path is possible for ACDP, but in the case of TMT, two different secondary reaction pathways are possible (vide supra).
Let us start with the reactant, S-c, which goes throughout one TS (S-TS co ) to the open diradical structure S-o (Figure 2).This S-TS co has a much lower Ea than the above S-TS X .The Ea for S-TS co of TMT turns out to be 17.7 and 18.3 kcal/mol for axial and equatorial methyl position isomers, respectively (Tables 3 and 4); and for ACDP Ea = 18.5 kcal/mol (Figure 2, Table 5).Considering the Ea as an average value (Ea aver ) of the axial and equatorial isomers, a function with the number (n) of methyl groups of the different derivatives of the series can be expressed by a linear equation Ea aver = 16.15 + 0.61n, R 2   Table 3.Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies (ν) b , spin contamination (S 2 ), and relative energies (∆E) c of critical points of the singlet-ground-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers.Acetaldehyde as the first intermediate product (Figure 3a) and acetone as the second product (Figure 4b).

Structure
Steps  From now on, one of the products of the step (either acetaldehyde or acetone), S-b′, is removed, and the oxide-peroxide diradical intermediate (S-b) continues reacting.A second molecule of acetaldehyde or acetone and one oxygen molecule will be the products of this reaction step (Tables 3-5).When S-b is •O-CH(CH3)-O-O•, TMT has the same Ea of 11.6 kcal/mol (as average value).This value is similar to MFDP, and DMT [8,9].The Er of and DMT, the rate-limiting step looks like the second step, where S-TSob′ has the highest energy barrier [7][8][9].However, these energy differences are very small and both steps can be considered limiting step at experimental conditions.Our experimental Ea for ACDP obtained in this work is 22.7 kcal/mol, which is consistent with our theoretical value of 18.5 kcal/mol of S-TSco (Table 5).4. Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies (ν) b , spin contamination (S 2 ), and relative energies (∆E) c of critical points of the singlet-ground-state PES' of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers.Acetone as intermediate product of step 2 (Figure 3b) and acetaldehyde as product of the last step (Figure 4a).
The next step is the TS of a broken C-O bond, S-TS ob ′ , which depends on either the twomethyl-or single-methyl-substituted C atom at TMT (Figure 3).In the first case, acetone will be obtained as the product (S-b ′ , Branch SII-ii of Scheme 1), while in the second case, the mono-substituted carbon (R = H, CH 3 in Scheme 1) will give acetaldehyde as the first product of reaction (S-b ′ , Branch SII-i of Scheme 1).Let us start with the mono-substituted carbon (Branch SII-i, Scheme 1).In Figure 3a the S-o intermediate goes throughout S-TS ob ′ to S-b ′ (second intermediate product).The Ea of S-TS ob ′ is 17.1 and 17.5 kcal/mol for axial and equatorial isomers with respect to the S-o intermediate.In S-TS ob ′ , the distinction between equatorial and axial isomers of CH 3 positions is not clear.Nonetheless, some slightly structural differences are found in the potential energy surface (PES) critical points.For this reason, average values are also taken.The transition vectors are mainly formed by the vibrations of C•••O distance-which are 1.704 and 1.707 Å for axial and equatorial isomers, respectively-and the spin density is shared by all oxygen atoms.The product is an acetaldehyde molecule plus an oxide-peroxide diradical, S-b ′ (Figure 3), whose Er is exothermic with respect to the S-o intermediate.However, it is slightly endothermic with respect to S-c.
When acetone is left as the first product in TMT (Branch ii in Scheme 1), the corresponding Ea of TS ob ′ and Er ob ′ are 15.4 and −13.9 kcal/mol, respectively, for the axial isomer and 15.4 and −12.7 kcal/mol, respectively, for the equatorial isomer (Figure 3b).This reaction is more favorable than the above, where the aldehyde is left out, though the differences are not meaningful.ACDP yields acetone and a peroxyacetone diradical as reaction products in the second stage, giving 15.31 and −14.83 kcal/mol for Ea and Er, respectively.
Considering the Ea as an average value (Ea aver ) of the axial and equatorial isomers and the two possible intermediates formed in this step of the different derivatives of the series, a linear function can be obtained.These Ea values low down to 1.
The slope is quite small, indicating that this TS is quite constant with respect to the number of methyl groups.The average values of axial and equatorial isomers give an exothermicity of 2.5 kcal/molCH 3 with respect to S-o, resulting from the linear equation Er ob ′ = −5.0− 2.5n (R 2 = 0.98).The correlation is low because the product's disposition is not very regular.In general, the Ea of this reaction step decreases with the number of methyl groups, and Er is growing in exothermicity.These facts indicate the additivity of methyl groups in our methyl-tetroxane derivatives; that is, each methyl group contributes to decreasing the activation energy.
From now on, one of the products of the step (either acetaldehyde or acetone), S-b ′ , is removed, and the oxide-peroxide diradical intermediate (S-b) continues reacting.A second molecule of acetaldehyde or acetone and one oxygen molecule will be the products of this reaction step (Tables 3-5).When S-b is •O-CH(CH 3 )-O-O•, TMT has the same Ea of 11.6 kcal/mol (as average value).This value is similar to MFDP, and DMT [8,9].The Er of the acetaldehyde product (Scheme 1, Branch SII-ii) is −29.01 and −28.64 kcal/mol with respect to b for TMTax and TMTeq, respectively (Figure 4a, Table 4).However, when the acetone is the second product, S-b is where the Ea and Er are 10.73 and −31.10 kcal/mol, respectively, for TMT and ACDP (Figure 4b).Therefore, the reaction of TMT and ACDP in the singlet state is exothermic in the last step of reaction.
In the derivatives with three and four methyl groups, the largest TS is that of the first step-that is, S-TS co , which could be the rate-limiting step.On the contrary, in FDP, MFDP, and DMT, the rate-limiting step looks like the second step, where S-TS ob ′ has the highest energy barrier [7][8][9].However, these energy differences are very small and both steps can be considered limiting step at experimental conditions.Our experimental Ea for ACDP obtained in this work is 22.7 kcal/mol, which is consistent with our theoretical value of 18.5 kcal/mol of S-TS co (Table 5).

Thermolysis Reaction as Triplet State (Scheme 2)
For ACDP, a quasi-chair excited structure (T-qc) in the PES of the triplet state has been found at 14.9 kcal/mol above S-c.This structure is not as symmetric as S-c and has a From T-o a transition state T-TS ob ′ is also found.This new structure is the TS for breaking and forming the first acetone/acetaldehyde molecule.In TMT, different values are found depending on if either acetone or acetaldehyde is the first product in b ′ .If acetaldehyde is the first product in b ′ (Scheme 2, Branch i), Ea is 16.98 and 15.13 kcal/mol for the axial and equatorial isomers, respectively (Figure 5a).If acetone is the first product (T-b ′ , Scheme 2, Branch ii), Ea is 15.20 kcal/mol for both axial and equatorial isomers (Figure 5b).In ACDP, Ea is 15.13 kcal/mol with respect to T-o (Figure 5b).The Ea as a function of the number of methyl groups from FDP to DMT was previously calculated with a quadratic function [9].With the addition of the three and four substituted derivatives to the system, the quadratic function is Ea = 25.3 − 6.2n + 0.9n 2 (R 2 = 0.97), similar to that found in Reference [9] (Tables 6-8).Table 6.Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies (ν) b , spin contamination (S 2 ), and relative energies (∆E) c of critical points of the triplet-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers.Acetaldehyde as the first intermediate product (Figure 5a) and acetone as the second product (Figure 6a).

Structure
Step

C•••O distances for T-TS ob
′ are 1.708 Å for ACDP and 1.704 and 1.715 Å for TMT as average values of the position axial/equatorial isomers when acetaldehyde and acetone are the first products, respectively.These distances are not quite different from those of S-TS ob ′ .TV (Figure 5) is viewed as a vibration of the C•••O breaking/forming the C•••O bond.These distances vary with the number (n) of methyl groups with a linear function: d C•••O = 1.714 − 0.0014n (R 2 = 0.991).The slope of the function is very small, where the structure of the TS can be considered very stable as the number of methyl groups increases.This equation is very close to the d C•••O equation of the S-TS ob ′ , which indicates an identity of structures in the TS ob ′ of the diradical part of reaction independently of the multiplicity of the state.Additionally, the o diradicals also have very close structures for both multiplicity states, as already mentioned above.From this, we can extract that if both multiplicity states have equal reactants and TS', the reaction path should be very close.The spin density of T-TS ob ′ is shared over all oxygen atoms in the structure (Figure 5).

Structure
Step  T-TS ob ′ goes to T-b ′ , a super-molecular system with an acetone/acetaldehyde molecule and a peroxy-acetaldehyde/acetone diradical, with a C acetone/acetaldehyde •••O peroxy distance of 3.956 and 3.473 Å for TMT aver(ax-eq) and acetone as the first product of the reaction; 3.647 Å is the average of both position isomers when the acetaldehyde is the first product of the reaction.The reaction energy is −11.06 kcal/mol with respect to T-o for TMT as an average value of both position isomers when acetaldehyde is the first product of reaction (Figure 5a, Table 6); and −13.36 kcal/mol for the equivalent value when acetone is the first product of the reaction (Figure 5b, Table 7).The Er as a function of the number of methyl groups behaves decreasing with the linear function Er = −5.5 − 2.9n, R 2 = 0.93.The spin density of the triplet state is shared on the three oxygen atoms of the peroxy-acetone radical.From now on, the acetone molecule is removed from the super-molecular system and the critical point is optimized from the peroxy-acetone radical.Nonetheless, the energy of acetone molecule is constantly added up to the energy of the critical points in order to preserve the whole energy of the system.Following this diradical intermediate, T-b, a T-TS bp yields the final product (T-p, Figure 6).The Ea is really the lowest of the T-PES and S-PES with values of 6.88-7.90kcal/mol (Tables 6-8).The Ea as a function of the number of methyl groups is also linear: Ea = 9.48 − 0.67n, R 2 = 0.98.In the T-PES all Eas low down with the increasing number of methyl groups, indicating the reaction rate is faster with the methyl substitution.Besides, this TS bp is the lowest TS of the S-and T-PES.The last TS gives rise to the second acetaldehyde/acetone molecule plus the 3 O 2 molecule.The TVs clearly show this formation (Figure 6a,b).The spin density also is shared over the three oxygen atoms of the system.The highest Ea on T-PES is that of T-TS ob ′ , being the rate-limiting step of the reaction, which agrees with the other derivatives of the series.From T-o a transition state T-TSob′ is also found.This new structure is the TS for breaking and forming the first acetone/acetaldehyde molecule.In TMT, different values are found depending on if either acetone or acetaldehyde is the first product in b′.If acetaldehyde is the first product in b′ (Scheme 2, Branch i), Ea is 16.98 and 15.13 kcal/mol for the axial and equatorial isomers, respectively (Figure 5a).If acetone is the first product (T-b′, Scheme 2, Branch ii), Ea is 15.20 kcal/mol for both axial and equatorial isomers (Figure 5b).In ACDP, Ea is 15.13 kcal/mol with respect to T-o (Figure 5b).The Ea as a function of the number of methyl groups from FDP to DMT was previously calculated with a quadratic function [9].With the addition of the three and four substituted derivatives to the Finally, the products are 3 O 2 plus the second acetone/acetaldehyde molecule.The reaction energy with respect to T-o is exothermic, with −58.2 kcal/mol for ACDP (Er ob ′ + Er bp , Table 8) and −54.2 and −53.9 kcal/mol for TMT aver(ax-eq) for acetone (Table 6) and acetaldehyde (Table 7) as the second product, respectively.A linear function is also found for this step of reaction, Er op = −41.4− 4.3n, R 2 = 0.998, describing the reaction energy of op, giving an exothermicity of 4.3 kcal/molCH 3 .Considering the first step of the reaction in singlet state for obtaining the biradical S-o, the reaction energies of the reactionin triplet state with respect to reactant S-c are −44.7 kcal/mol for ACDP and −41.4 and −40.7 kcal/mol for TMT aver(ax-eq) for acetone and acetaldehyde as final products, respectively.If the reaction would proceed in singlet state PES, the Er of the equivalent critical point for S-p, with respect to S-c is −24.0 kcal/mol and −29.9 kcal/mol for TMT aver(ax-eq) for acetone and acetaldehyde as final products, respectively and −32.4 kcal/mol with respect to S-c for ACDP (Tables 3-5).This indicates that the reaction energies on T-PES are much more exothermic than in the S-PES, as expected based on the fact that the singlet molecular oxygen is more unstable than triplet molecular oxygen.The results in the T-PES are more in agreement with the experimental reaction energies of the diperoxide compounds.
The influence of the methyl groups in the thermolysis reaction of tetroxane derivatives increases the exothermicity of the reaction, and they change the rate-limiting step from the second step to the first step, although the differences in Ea are very small, and both TSs could be viewed as limiting steps at experimental conditions; additionally, they also change the geometry of intermediates and transition states in an approximately linear way.The singlet and triplet states present very close structures at the different critical point intermediates of the reaction.These similarities of structures indicate a possible nonadiabatic transition from the singlet to triplet states in such a way that reaction would be more exothermic if this transition of different multiplicity states is produced, which would agree with the high exothermicity of the tetroxanes.Therefore, it is worth providing general details on the intersystem crossing process, as described in the next section.

Excited States and Non-Adiabatic Chemistry of Tetroxane and the Methylated Systems
The lowest-lying singlet and triplet excited states of tetroxane (FDP) at the reactant structure correspond to electronic excitations from the lone pairs of the oxygen atoms (n i ) or the σ bonding orbitals of the -O-O-bonds (σ i ) to the σ anti-bonding orbitals of those bonds (σ * j )-that is, n i → σ * j or σ i → σ * j excitations (see orbitals in Figure 7).These electronic configurations are associated with the weakening of the peroxide bonds.The vertical energy position of the three first singlet (S i ) excited states (all of them n i → σ * j in nature) relative to the energy of the ground state (S 0 ) ranges from 5.09 to 5.78 eV (Table 9).The manifold of triplet states (also n i → σ * j in nature) is lower-lying than that of the singlet states, being in a range of 3.87-4.81eV for the first four of them.Test computations show that σ i → σ * j configurations appear on the geometry of the reactant at higher energies (Table 9).Upon -O-O-bond stretching, the mentioned electronic configurations become stabilized as is confirmed by the calculated energy of the excited S i and T i states at the diradical geometry.Here, the four singlets and four triplet states appear in an energy range lower than 1 eV (Figure 8).T 1 has the same energy as the ground state (S 0 ), and other states (S 1 , S 2 , T 2 , and T 3 ) are located only 0.3 eV above S 0 .Considering that the kinetic energy released after the transition state of the -O-O-bond breaking is of the same order; those excited states may be populated in the non-adiabatic dynamical process, which follows the peroxide bond decomposition.To analyze the efficiency of the non-adiabatic transitions at the diradical region, SOCs were calculated between the different singlet and triple states (Table 10).The coupling between S0 and T1 is negligible, and the maximum value is obtained with T3, which is located at 0.29 eV.Additionally, a spin orbital coupling (SOC) between S2 and T1 at 0.29 eV also has a very high value.Both SOCs indicate two possibilities for transferring the To analyze the efficiency of the non-adiabatic transitions at the diradical region, SOCs were calculated between the different singlet and triple states (Table 10).The coupling between S 0 and T 1 is negligible, and the maximum value is obtained with T 3 , which is located at 0.29 eV.Additionally, a spin orbital coupling (SOC) between S 2 and T 1 at 0.29 eV also has a very high value.Both SOCs indicate two possibilities for transferring the population from S 0 to T 1 along the thermal decomposition: (i) S 0 transfers the population initially to T 3 via strong SOC and subsequently by an ultrafast internal conversion process T 3 decays to T 1 , or (ii) S 0 transfers the population initially to S 2 via an internal conversion process and next the molecule changes to the T 1 manifold through strong SOC between S 2 and T 1 .Analysis of the open-shell spin densities at the diradical region shall help to provide a rationale on the distinct efficiencies related to the population transfer between the singlet and triplet states (Figure 9).Such densities reflect the distinct occupation of the unpaired electrons of the diradical in the oxygen 2p x and 2p y atomic-like orbitals, which are perpendicular to the C-O bonds.As can be seen, high SOC values are obtained when singlets and triplets have one unpaired electron in a perpendicular 2p orbital, in accordance with El-Sayed's rules [18].
Once we comprehend the electronic-structure features of the parent tetraoxane molecule (FDP), it is worth extending the analysis to the methylated molecules of the family, Figure 8 shows that ACDP has the electronic states closer in energy compared with FDP.Energies of S 3 and T 4 show a clear trend of decreasing energy upon methylation.Regarding the most accessible states (S 1, S 2 , S 3 , T 2 , and T 3 ), we see that an odd number of methyl group substitutions (MFDP and TMT) increases the energy gap, while those states are energetically slightly more accessible for FDP and the molecules with an even number of methyl groups.Comparison of SOCs in the distinct molecules of the family (Table 10) also do not show an efficient transition from S 0 to T 1 at the diradical region and a higher efficiency for S 0 -> T 3 and S 2 -> T 1 .The highest values appear also for the molecules with an even methylation.ACDP, which shows the lowest energy position of the S 2 and T 3 states at the diradical region and high SOCs arises as the most efficient system in the family for non-adiabatic decomposition.The dynamical effects, not considered in this work, are expected to favor this triplet population even more, as inferred from the study by Vacher et al. comparing the dynamics of 1,2-dioxetane and the tetramethylated compound [19].By increasing the mass of the substituents, the time spend at the diradical region increases, allowing a more effective re-distribution of the population among the near-degenerate singlet and triplet states and accordingly higher yields for the evolution on T 1 [19].These transitions are more favorable for FDP and even-methylated derivatives (the highest SOCs and the least ∆E), being less efficient at the odd-methylated compounds.For FDP and the even-methylated derivatives, SOC is linear with the number of methyl groups (SOC S2->T1 = 96.15-0.96•n;R 2 = 0.999; SOC of DMT being the average of the three position isomers); SOC of S 0 to T 3 equally decreases with the number of methyl groups, but non-linearly.∆E (Figure 8) of the tetramethyl derivative is the least of all derivatives.SOCs of both transitions of the two odd-methylated derivatives increase with the number of methyl groups, SOCs at S 2 to T 1 being the highest.Analysis of the open-shell spin densities at the diradical region shall help to provide a rationale on the distinct efficiencies related to the population transfer between the singlet and triplet states (Figure 9).Such densities reflect the distinct occupation of the unpaired electrons of the diradical in the oxygen 2px and 2py atomic-like orbitals, which are perpendicular to the C-O bonds.As can be seen, high SOC values are obtained when singlets and triplets have one unpaired electron in a perpendicular 2p orbital, in accordance with El-Sayed's rules [18].Once we comprehend the electronic-structure features of the parent tetraoxane molecule (FDP), it is worth extending the analysis to the methylated molecules of the family, Figure 8 shows that ACDP has the electronic states closer in energy compared with FDP.Energies of S3 and T4 show a clear trend of decreasing energy upon methylation.Regarding the most accessible states (S1, S2, S3, T2, and T3), we see that an odd number of methyl group substitutions (MFDP and TMT) increases the energy gap, while those states are energetically slightly more accessible for FDP and the molecules with an even number of methyl groups.Comparison of SOCs in the distinct molecules of the family (Table 10) also do not show an efficient transition from S0 to T1 at the diradical region and a higher efficiency for S0 -> T3 and S2 -> T1.The highest values appear also for the molecules with an even methylation.ACDP, which shows the lowest energy position of the S2 and T3 states at the diradical region and high SOCs arises as the most efficient system in the family for

Materials
DMT was prepared by reacting 0.6 mL of acetaldehyde, which is added to a stirred and cooled solution at 263.15 K containing 2.5 mL of ethanol (absolute; Merck, Darmstadt, Germany), 0.3 mL of hydrogen peroxide (56% v/v) and 1.5 mL of sulfuric acid (98%, Merck).The mixture was stirred, and sulfuric acid was added up to 3 mL, precipitating a white solid.Then, water was slowly added to complete a volume of 5 mL in order to increase the insolubility of the peroxide precipitate in the alcoholic medium.A white precipitate with a gelatinous characteristic was obtained, which was washed and centrifugated repeatedly with distilled water until neutral pH.Thus, the wet solid was dissolved directly in benzene (avoiding the incorporation into the organic phase of possible polar impurities).This dissolution was analyzed using the gas chromatographic (GC) method by using an Agilent 7890-A chromatograph (Santa Clara, CA, USA).
ACDP was synthesized by carefully adding acetone to a stirring solution of hydrogen peroxide in strong acid medium at 253 K, as it was previously proposed [7].ACDP solid was characterized by IR spectroscopy by using a Nicolet spectrometer (Thermo Fisher Scientific, Waltham, MA, USA), and a solution was analyzed by UV spectroscopy with a CampSpec M330 spectrophotometer (Leeds, UK) and GC.THF was used as a solvent in the thermolysis reaction of ACDP.

Analytic and Kinetic Methods
The GC method was used to study the kinetic of the thermolysis reaction of DMT and ACDP.This technique is very suitable for small quantities of samples.The reactor in this method is the injection chamber of a gas chromatograph, where the flow and reaction time (40 s) are maintained constant, and the reaction temperature is a parameter.When the compound flows out of the reactor and goes to the chromatographic column, temperature lowers enough to stop the thermolysis reaction, but temperature is maintained high enough to produce the chromatographic separation [12] (∆T = −150 K in the most adverse case).The method was also used for FDP, obtaining good results [7].
A solution (1 µL) of DMT (0.001 M) in benzene was introduced in the injection chamber.The reaction was controlled by the injector temperature within the range of 493.15-543.15K, and the rate of ∆T = 10 K.As the flow and linear rate were kept constant, the reaction time is consequently kept also constant.The reaction products and reactant were analyzed by GC with an electronic device of constant flow (Q = 1.40 mL/min; V L = 20 cm/s) using nitrogen as carrier gas.A flame-ionization detector and a melted silica capillary column (25 m × 0.53 mm) coated with 5% of phenyl-methylpolysiloxane were used.Temperature of the chromatographic column was started at 313 K for 2 min, then programmed at a rate of 30 K/min up to 463 K, and held constant for 10 min.
A solution (1 µL) of ACDP (0.02 M) in THF was introduced in the injection chamber under the same experimental conditions as DMT.The reaction was controlled by the injector temperature.
Residence time, which coincides with the reaction time, in the injector chamber was calculated by Reference [12]: t = l/v (5) where l (cm) is the liner length and v is the linear rate of the carrier gas.All previous works have found the first-order reaction model for this kind of thermolysis reaction [7,12].Then, Ea was calculated by the Arrhenius equation in a first-order integral kinetic equation such as: ln[ln(C 0 /C)] = −Ea/RT + ln(A•t) (6) where C 0 and C are the initial and remaining concentrations of the reactant, respectively; T is the temperature; A is the frequency factor of the Arrhenius equation; and t is the residence time, Equation (4).ln[ln(C 0 /C)] is fitted as a function of 1/T to a linear equation.

Computational Methods
The electron density properties of the studied molecules (ACDP and TMT) were calculated according to references [7][8][9].The geometry and energy of the reactant, intermediates, and products were computed by using the Gaussian 09 program [20], at the DFT/BHANDHLYP/6-311 + G** level [21][22][23].Critical points (minimums and TS) of the S-and T-PES were calculated independently and optimized with the Berny method [24].The topology of the critical point was characterized by means of the harmonic frequencies.Transition vectors (TVs) for the TSs were inspected and drawn with the help of the GaussView package [25].The intrinsic reaction coordinate method was used in order to connect the reaction paths between the TS and reactant, products or intermediates [26,27].The relative energies of the critical points of the PES were calculated by adding the total energy to the zero-point energy (ZPE).
The tetroxane series (FDP, MFDP, DMT, TMT, and ACDP) with their relevant axial and equatorial conformers, as determined in the current and previous related works, was revisited to obtain a high-level characterization of distinct nonadiabatic transitions between T and S states occurring along the thermal decomposition reaction.Four S states (S 0 -S 3 ) and four T states (T 1 -T 4 ) were energetically characterized accounting for the spin orbit coupling

Scheme 1 .
Scheme 1. Reaction paths of the singlet ground state ACDP and TMT by single-step/concerted (Branch S-I: from S-c to S-P) and stepwise (Branch S-II: from S-c to S-p) mechanisms.Arrow up and down means spin up and down, respectively.

Scheme 1 .
Scheme 1. Reaction paths of the singlet ground state ACDP and TMT by single-step/concerted (Branch S-I: from S-c to S-P) and stepwise (Branch S-II: from S-c to S-p) mechanisms.Arrow up and down means spin up and down, respectively.

Figure 1 .
Figure 1.(a) Concerted mechanism of S-PES for TMT and ACDP, black color for ACDP, blue for axial TMT, and green for equatorial TMT, and blue arrows indicate the transition vector of TS.(b) Activation energy of TS of the concerted mechanism (TSX) as a function of the number of methyl groups in FDP.Values for FDP, MFDP, and DMT come from references [7-9].

Figure 1 .
Figure 1.(a) Concerted mechanism of S-PES for TMT and ACDP, black color for ACDP, blue for axial TMT, and green for equatorial TMT, and blue arrows indicate the transition vector of TS.(b) Activation energy of TS of the concerted mechanism (TS X ) as a function of the number of methyl groups in FDP.Values for FDP, MFDP, and DMT come from references [7-9].
= 0.993.The most important distance in the TS, and so the simplified reaction coordinate, in the structure of S-TS co , is the O•••O distance, which, in the average of the axial and equatorial isomers (d O•••Oaver ), also shows a linear function with the number of methyl groups: d O•••Oaver = 1.90 + 0.014n, R 2 = 0.96.Indeed, the Ea and d O•••O as a function of the number methyl groups show qualitative parallel trends, increasing the distance as increasing Ea, but with a much lower slope of 0.014 Å/CH 3 , indicating that the main coordinate of the TS has a small variation on the effect of substituents in the same reaction.state.This intermediate gives endothermic reaction energy, increasing with the number (n) of methyl groups, following also a linear function: Erco = 9.0 + 1.22n; R 2 = 0.992.The slope of Erco is twice the slope of Ea-TSco, indicating that the methyl substitution effect is twice as large in the diradical open intermediates as in TSco.Indeed, the linear equation relating both Erco and Ea (TSco) is Erco = −23.5 + 2.0Ea (R 2 = 0.996).The O•••O distance of S-o follows the linear equation dO•••O = 3.519 − 0.074n, R 2 = 0.997.Additionally, the reaction coordinate O•••O of S-o intermediate decreases with a slope of 4.6 Å with respect to this O•••O distance in S-TSco, following the linear equation dO•••Oo = 12.2 − 4.6d(O•••O)TSco (R 2 = 0.99).All these effects agree with the Leffler-Hammond postulate[16,17].

Figure 2 .
Figure 2. First step of stepwise mechanism S-PES of TMT and ACDP.Black color for ACDP, blue for axial TMT, and green for equatorial TMT.The main non-bonding distances are shown in Å.The spin density is indicated in brackets.Blue arrows indicate the transition vector of TS.

Figure 3 .
Figure 3. Second step of the stepwise mechanism: (a) when acetaldehyde is the first product for both axial and equatorial isomers of TMT; (b) for ACDP and for TMT when acetone is the first product for axial and equatorial isomers.The interatomic distance of the breaking bonds are shown, and the spin densities are in brackets.The values of ACDP are in black, and those of TMT are in blue for axial TMT and green for equatorial TMT.Blue arrows mean the transition vector of TS.

Figure 3 .
Figure 3. Second step of the stepwise mechanism: (a) when acetaldehyde is the first product for both axial and equatorial isomers of TMT; (b) for ACDP and for TMT when acetone is the first product for axial and equatorial isomers.The interatomic distance of the breaking bonds are shown, and the spin densities are in brackets.The values of ACDP are in black, and those of TMT are in blue for axial TMT and green for equatorial TMT.Blue arrows mean the transition vector of TS.

Figure 4 .
Figure 4. Third step between S-b to S-p of the S-PES: (a) for TMT for axial and equatorial isomers when acetaldehyde and oxygen are the last products; (b) for ACDP and for TMT for axial and equatorial isomers when acetone and oxygen are the last products.The interatomic distance of the breaking bond is indicated in Å, and the electron densities are shown in brackets.Black color for ACDP, blue for axial TMT, and green for equatorial TMT.Blue arrows mean the transition vector of TS.Table3.Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies (ν) b , spin contamination (S 2 ), and relative energies (ΔE) c of critical points of the singlet-ground-state PES of the thermolysis reaction of 3,3,6-trimethyl-1,2,4,5-tetroxane (TMT) in position C-6 axial (ax) and equatorial (eq) isomers.Acetaldehyde as the first intermediate product (Figures 3a) and acetone as the second product (Figure4b).

Figure 4 .
Figure 4. Third step between S-b to S-p of the S-PES: (a) for TMT for axial and equatorial isomers when acetaldehyde and oxygen are the last products; (b) for ACDP and for TMT for axial and equatorial isomers when acetone and oxygen are the last products.The interatomic distance of the breaking bond is indicated in Å, and the electron densities are shown in brackets.Black color for ACDP, blue for axial TMT, and green for equatorial TMT.Blue arrows mean the transition vector of TS.

Table
The first intermediate in the reaction is the diradical open structure, S-o, with a distance between both O atoms of 3.22 Å and a spin density of −0.94 and 0.94 e − (Figure2), indicating that both electrons have different spins, maintaining the identity of the singlet state.This intermediate gives endothermic reaction energy, increasing with the number (n) of methyl groups, following also a linear function: Er co = 9.0 + 1.22n; R 2 = 0.992.The slope of Er co is twice the slope of Ea-TS co , indicating that the methyl substitution effect is twice as large in the diradical open intermediates as in TS co .Indeed, the linear equation relating both Er co and Ea (TS co ) is Er co = −23.5 + 2.0Ea (R 2 = 0.996).The O•••O distance of S-o follows the linear equation d O•••O = 3.519 − 0.074n, R 2 = 0.997.Additionally, the reaction coordinate O•••O of S-o intermediate decreases with a slope of 4.6 Å with respect to this O•••O distance in S-TSco, following the linear equation d O•••Oo = 12.2 − 4.6d(O•••O) TSco (R 2 = 0.99).All these effects agree with the Leffler-Hammond postulate 3 kcal/molCH 3 , following the linear equation Ea(TS o b ′ ) = 20.6 − 1.3n (R 2 = 0.93).The C•••O distances have a linear function with respect to the number (n) of methyl groups: 82 Å, much larger than the other peroxide bonds (1.41 Å).The methyl groups are in the same conformation as S-c.This excited structure has a spin density of 0.942 e − on the two open oxygen atoms.The equivalent structure for TMT has not been found.The second structure found in the T-PES is the open diradical structure (T-o) close to S-o.This T-o presents an energy of 13.7 kcal/mol with respect to S-c and only −1.2 kcal/mol with respect to T-qc.Then, the energy of T-o is very close to T-qc, where the quasi-open structure of T-qc has a large d O•••O distance similar to the open structure T-o.If a transition state exists from T-qc to T-o, we could not find it, but it should be with low Ea from T-o to T-qc.The energy and structure of T-o are very close to the S-o, having a spin density of 0.939 e − on each open oxygen of the peroxide bond, and an O•••O distance of 3.222 Å for ACDP, and 3.286 and 3.330 Å for TMT axial and equatorial isomers, respectively (Figure 5).The O•••O distances decrease as the number (n) of methyl groups following the linear equation d O•••O = 3.527 − 0.075n (R 2 = 0.998).This equation and distances are very close to the S-o-d O•••O distances, indicating the similar o structures for the different derivatives and their very close variation with the number of methyl groups.

Figure 5 .
Figure 5.The first steps of the stepwise mechanism between T-o and T-b′ of the T-PES: (a) TMT when acetaldehyde is the first product for both axial (ax) and equatorial (eq) isomers; (b) TMT and ACDP when the acetone is the first product.Black color is for ACDP, blue for axial TMT, and green for equatorial TMT.Blue arrows indicate the transition vector of TS, the breaking bond distances are indicated in Å, and the spin density values are in brackets.

Figure 5 .
Figure 5.The first steps of the stepwise mechanism between T-o and T-b ′ of the T-PES: (a) TMT when acetaldehyde is the first product for both axial (ax) and equatorial (eq) isomers; (b) TMT and ACDP when the acetone is the first product.Black color is for ACDP, blue for axial TMT, and green for equatorial TMT.Blue arrows indicate the transition vector of TS, the breaking bond distances are indicated in Å, and the spin density values are in brackets.

Figure 6 . 6 .
Figure 6.The second step between T-b to T-p of the T-PES of (a) TMT and ACDP when the acetone and oxygen are the products; (b) TMT when the acetaldehyde and oxygen are the products.The Figure 6.The second step between T-b to T-p of the T-PES of (a) TMT and ACDP when the acetone and oxygen are the products; (b) TMT when the acetaldehyde and oxygen are the products.The distance of the breaking bond is indicated in Å, and the spin density values are in brackets.The TVs are indicated by blue arrows.

Molecules 2024 ,Scheme 2 .
Scheme 2. Reaction paths of the triplet ground states of ACDP and TMT.The two arrows down means two parallel spins.

Scheme 2 .
Scheme 2. Reaction paths of the triplet ground states of ACDP and TMT.The two arrows down means two parallel spins.

Molecules 2024 , 27 Figure 8 .
Figure 8. Vertical energies of the four singlet and four triplet electronic states at the ground-state geometry of the diradical obtained with the CASSCF/CASPT2(12-in-8)/ANO-L-VDZP calculations.Green levels are excited triplet states (from T1 to T4), and blue levels are ground singlet state (S0) and their excited states (S1 to S3) of each derivative in increasing order of energy.

Figure 8 .
Figure 8. Vertical energies of the four singlet and four triplet electronic states at the ground-state geometry of the diradical obtained with the CASSCF/CASPT2(12-in-8)/ANO-L-VDZP calculations.Green levels are excited triplet states (from T 1 to T 4 ), and blue levels are ground singlet state (S 0 ) and their excited states (S 1 to S 3 ) of each derivative in increasing order of energy.

Figure 9 .
Figure 9. CASSCF (12-in-8)/ANO-L-VDZP electron spin density for the triplet states of the diradical of FDP.The corresponding densities for the singlets S 0 , S 1 , S 2 , and S 3 are similar to those for T 4 , T 2 , T 3 , and T 1 , respectively.

Table 1 .
Temperature (K) and rate constant (k) (s −1 ) of DMT 0.001 M in solution of benzene and ACDP 0.02 M in solution of tetrahydrofuran (THF).

Table 2 .
Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies ν b , spin contamination S 2 , and relative energies (ΔE) c of critical points of the singlet-ground-state PES of the thermolysis reaction of TMT and ACDP.Concerted mechanism.

Table 2 .
Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies ν b , spin contamination S 2 , and relative energies (∆E) c of critical points of the singlet-ground-state PES of the thermolysis reaction of TMT and ACDP.Concerted mechanism.Imaginary frequencies in cm −1 .c Relative energies in kcal/mol, ∆E = [E total + ZPE](c) − [E total + ZPE](i), where i, represents any critical point of the species in the table.d In position C-6 axial (ax) and equatorial (eq) isomers.
a Hartree.b
Imaginary frequencies in cm−1 .cRelative energies in kcal/mol, ∆E = [E total + ZPE](o/b) − [E total + ZPE](i), where i represents any critical point of the species in the table.d The relative energies with respect to Sc are 28.78 and 31.03kcal/mol for the axial and equatorial isomers, respectively.
a Hartrees.b e The relative energies with respect to Sc are −3.67 and −4.33 kcal/mol for the axial and equatorial isomers, respectively.
represents any critical point of the species in the table.dTherelative energies with respect to Sc are 27.16 and 28.95 kcal/mol for the axial and equatorial isomers, respectively.e The relative energies with respect to Sc are −2.16 and −2.12 kcal/mol for the axial and equatorial isomers, respectively.

Table 8 .
Total energies (E) a , zero-point energies (ZPE) a , imaginary frequencies (ν) b , spin contamination (S 2 ), and relative energies (∆E) c of the critical points of the triplet PES for the thermolysis reaction of 3,3,6,6-tetramethyl-1,2,4,5-tetroxane (ACDP).a Hartrees.b Imaginary frequencies in cm −1 .c Relative energies in kcal/mol, ∆E = [E total + ZPE](o/b) − [E total + ZPE](i), where i represents any critical point of the species in the table.d With respect to T-o. e The relative energy with respect to Sc is 28.70 kcal/mol.f The relative energy with respect to Sc is −1.56 kcal/mol.

Table 9 .
Spectroscopic characterization of the reactant.E VA corresponds to the CASPT2 (12-in-8)/ANO-L-VDZP vertical absorption energy (eV) of FDP.Main nature of the transitions and the corresponding weight of the electronic configurations are also shown (see orbitals in Figure7).

Table 10 .
SOC complex vector norms (see Computational Details) in cm −1 for the tetraoxane series.Data obtained with the CASSCF(12,8)/ANO-L-VDZP method.Note that T i encompasses the three m s components of the triplet state.